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Reference |1| (henceforth referred to as /) considered a general class of problems involving the evolution of large 
systems of globally coupled phase oscillators. It was shown there that, in an appropriate sense, the solutions to these 
problems are time asymptotically attracted toward a reduced manifold of system states (denoted M). This result has 
considerable utility in the analysis of these systems, as has been amply demonstrated in recent papers [2|-[l3|- In this 
note, we show that the analysis of / can be modified in a simple way that establishes significant extensions of the 
range of validity of our previous result. In particular, we generalize / in the following ways: (1) attraction to M is 
now shown for a very general class of oscillator frequency distribution functions (?(w), and (2) a previous restriction 
^I^ on the allowed class of initial conditions is now substantially relaxed. In particular, with respect to point (1), we now 

show that the result of / (derived there for the special case of Lorentzian g(a;)) actually applies for a very general 
distributions g{uj) described in the third paragraph below. 

To proceed, we recall the general class of problems described by Eqs.(l) to (4) of /, where the time evolution of 

the quantity H[t) [appearing in Eq.(4) of /] is determined from the time evolution of the order parameter r(t) [given 

^vj ' by Eq.(2) of /], and, as discussed in /, this determination of H from r can arise in various ways depending on the 

\ particular problem under consideration. The oscillator distribution function F{9,uj,t) is then expressed [Eqs.(5) and 

I""!' (6) of /] in terms of a decomposition involving functions F^{9,ui,t) and F_(9,ui,t), where the analytic continuation 

^^ \ of F+ (F-) into Im{6) > {Im{6) < 0) has no singularities and decays to zero as Im{9) -^ +oo {Im{6) — > — oo). F+ 

\^ is further decomposed [Eq.(8) of /] into two parts, F+ = F^ + F^, where F^ lies on the reduced manifold M and has 

P^ . dynamics described by Eq.(9) of /, while F+ satisfies 
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F+ = 0. (1) 
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JZ^ , As argued in /, the time asymptotic evolution of the order parameter will be completely described by dynamics on 
M (i.e., by F^) provided that 

cn 
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in' lim / F+(e,uj,t)g(uj)duj = 0, (2) 

where g{ui) is the oscillator frequency distribution function. In / it is shown that Eq.® is satisfied for the case 
where gito) is Lorentzian and the analytic continuation of the initial condition, F^{9, io,0), into Im{uj) < has no 
singularities and approaches zero as |a;| -^ cx), Im{u!) < 0. Here we will show that ([5]) is satisfied much more generally. 

We consider a general class of oscillator frequency distributions (?(w) that are analytic for real ui and can be 
5-H I analytically continued into a strip S defined by > Im{uj) > ~a, cr > 0, in which g{uj) has no singularities and decays 
■ - ■ to zero as Re{uj) -^ ±oo fast enough so that J_ g{^r + iu)i)dLOr is well-defined. Some examples of g{u)) satisfying 
these conditions are the Maxwellian, g{uj) ^ exp{ — (a; — w)^/[2(Aa;)^]}, the Lorentzian, g{uj) ~ [(w — w)^ -|-(Aa;)^)]~"'^, 
and the sech distribution, g{ijj) ~ sech[(a; — Q)/ /S.lj\. In contrast, paper / restricted consideration to Lorentzian g{uj). 
We further assume that the initial condition F+(6',a;,0) satisfies the same condition in S as g{Lo). Note that this 
condition on F-^-{0, w, 0) is considerably weaker than the condition used in / where it was required that F^{6, u, 0) be 
analytic everywhere in Im{uj) < and F^{9, lj, 0) — >■ for |a;| — >■ oo, Im{uj) < 0. 

In order to show ([2]) for the above conditions on g and F+, we consider the solution of Eq.([T]) for a; in the strip S. 
Noting that w only appears as a parameter in ([1} we can, for the moment, regard uj — ujr + ioJi [uji < 0) as fixed. 
Furthermore, by means of the replacement, F+ -^ F+e"'"''*, we can, without loss of generality, take w^ — 0. Thus 
Eq.dl]) takes the same form as Eqs.(17) and (18) of /. Furthermore, it was shown [see Eqs.(19) to (31) of / and the 
accompanying discussion] that the solution of (18) and (19) of / approaches zero at i -> +oo. Thus, for any w G 5 
and Im{uj) < 0, we have that |F+| — > as i — >■ -|-cxd. To complete our argument, we now note that our conditions 
allow us to shift the integration path for the integral in ^ from the real w-axis to the line ujr + iuJi, > tUi > —a, 
— oo <ujr< +00. Since at every point w^ on this line [g(tij)| is bounded (because g{uj) has no singularities in S) and 
F+ — > as i — > oo, we conclude that Eq.([5]) is satisfied. 



The above establishes a greatly expanded range of situation in which Eq. (9) of / can be used for discovering the 
time asymptotic dynamics of phase oscillator systems. We note, however, that, while (9) of J is substantially simpler 
than the original problem, it is still nontrivial. Thus for the solution of (9) for the order parameter dynamics it may 
still be useful to employ special choices for g(w) such as the Lorentzian (see Ref.|3|). 

Finally, we note that our requirement that g{uj) be analytic means that our proof of attraction to M does not apply 
if g{io) is a delta function. Thus attraction to Af necessitates that g{uj) have a finite width. This is in line with the 
intuition that relaxation to M results from phase mixing of oscillators in a population with heterogeneous frequencies. 
The inapplicability of our result in the case where oscillators all have the same frequency is also consistent with past 
results [3, [ij,[ia| which imply that, when there is no frequency spread, the dynamics is not attracted to M. 
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